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An isolated classical chaotic system, when driven by the slow change of several parameters, re- 
sponds with two reaction forces: geometric friction and geometric magnetism. By using the theory 
of quantum fluctuation relations we show that this holds true also for open quantum systems, and 
provide explicit expressions for those forces in this case. This extends the concept of Berry curva- 
ture to the realm of open quantum systems. We illustrate our findings by calculating the geometric 
magnetism of a damped charged quantum harmonic oscillator transported along a path in physical 
space in presence of a magnetic field and a thermal environment. We find that in this case the 
geometric magnetism is unaffected by the presence of the heat bath. 

PACS numbers: 05.40.-a 03.65.Vf 



I. INTRODUCTION 

According to the adiabatic theorem, when a thermally 
isolated quantum system is subjected to an adiabatic 
cyclical driving it returns to the same state where it 
started from [l|. The only effect of the excursion is that 
the state acquires a phase. In 1984, Berry clarified that 
such a phase consists of a gauge-dependent (hence non- 
physical) part and a gauge-independent (hence physical 
and measurable) part, which is determined only by the 
geometry of the path described by the driving parame- 
ters @. This is by now customarily called Berry's phase. 
Soon later Hannay [|| and Berry [j] found the analogous 
classical phenomenon: when a classical intcgrable sys- 
tem undergoes an adiabatic cyclic evolution, the action 
variable remains constant, but the angle variable expe- 
riences an anholonomy effect, that is, it does not return 
to its original value but accumulates a shift, known as 
Hannay angle. More recently, Robbins and Berry [f| ad- 
dressed the question of whether there exists a geometric 
phase also in the case of classical chaotic systems. Their 
approach was based on the observation that the Berry 
phase is given by the flux of a two-form (Berry's curva- 
ture) through a surface bounded by the cyclic path in 
the parameter space. Therefore, they investigated the 
classical limit of the quantum two-form and found the 
expression 
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where Vift is the gradient in parameter space of the 
Hamiltonian evaluated at time t (i.e., the instantaneous 
negative force exerted by the external driving), x denotes 
cross product, (■ ■ - )e denotes microcanonical average at 
energy E and uj(E) is the corresponding density of states. 
Later Jarzynski Q showed that the surface integral of 
the classical chaotic two-form ([T]) measures a shift ac- 
cumulated along the chaotic trajectory on the constant 
energy hyper-surface, which generalizes the concept of 
Hannay's angle. In a subsequent paper Berry and Rob- 
bins Q re-derive Eq. ([1]) adopting a statistical mechan- 
ical approach. They considered an initially microcanon- 



ically distributed ensemble and focused on the average 
force with which the system reacts to the external driv- 
ing, F = — (Vif) E-, up to first order in the driving speed. 
They showed that the reaction force contains two terms: 
a friction-like force and a Lorentz-like force. The latter 
stems from a magnetic-like field, the so called geomet- 
ric magnetism, which is nothing but the classical chaotic 
two- form detailed in Eq. ([T]). The approach developed 
in Q is similar to Kubo's Linear Response Theory 0]. 
The main differences are: (i) in Kubo's theory the initial 
state is canonical whereas in Berry and Robbins theory it 
is microcanonical; (ii) Kubo's theory gives the response 
of a driven system to a weak perturbation, and gives, 
accordingly, a linear relationship between the response 
and the strength of the driving. The theory of Berry and 
Robbins instead yields the response of a driven system 
to a slow perturbation, and gives, accordingly, a linear 
relationship between the response and the speed of the 
driving. We might label Berry and Robbins theory Q a 
Microcanonical Adiabatic Linear Response Theory. 

In the recent years yet further attempts have been de- 
voted to investigate possible generalizations of Berry's 
phase. It has been realized that geometric phases can 
be used as effective and reliable tools for quantum com- 
putation [j| Thus, the study of geometric phases in 
realistic open quantum systems has become of paramount 
importance. This problem has been typically addressed 
from a dynamics point of view, that is, one researches 
proper definitions of phases with respect to non-unitary 
dynamics [TlTj2l1 | (which are the relevant ones for open 
systems) , instead of the unitary dynamics originally con- 
sidered by Berry. 

To the best of our knowledge, however, the problem 
of finding the geometric phase of an open quantum sys- 
tems was not previously addressed adopting a statistical 
mechanical method akin to the one employed by Berry 
and Robbins [3] ■ Here we pursue this task. Below we de- 
velop a Canonical Adiabatic Linear Response Theory for 
open quantum systems. As a main result we obtain the 
general expression of the field of geometric magnetism of 
open quantum systems, see Eq. (f21~j) below. The geo- 
metric phase is given by its surface integral in analogy 
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with the standard case. In developing the theory we take 
full advantage from the theory of quantum work fluctu- 
ation relations [13, HH , which can be formulated within 
two complementary viewpoints: exclusive and inclusive 
[24|,[25|. While the former is best suited to derive Kubo's 
Linear Response Theory (2(| the latter, as we will 
see below, is best suited to derive the searched Canonical 
Adiabatic Linear Response Theory. 



where AF = —ft' 1 1ii(Zr t /Zr ) is the difference of free 
energy of the total system states and g^ o . For O = 1 
Eg. ([7]) reproduces the quantum Jarzynski Equality [H, 
1301 ] - Eq. (J7J rnay be obtained from the nonequilibrium 
generating functional of Andrieux and Gaspard [27] by 
means of functional differentiation. Note that the free 
energy difference can be written in the following form 

& AF = -f^dR-(Q)% = -f T dm t -(Q)2. Using 
the notations 



II. ADIABATIC LINEAR RESPONSE OF OPEN 
QUANTUM SYSTEMS 

We consider a driven open quantum system in contact 
with a thermal bath at fixed temperature T = l/(fcg/3). 
Following the established procedure [28|, [2|| we hence 
close the system by coupling it to a thermal bath. The 
total (system+bath) Hamiltonian reads 



H(-R t ) = H B +H SB + H(R t ), 



(2) 



where Hb is the bath Hamiltonian, H$b is the (possibly 
strong) system/bath coupling, and if (Rt) is the system 
Hamiltonian: 



H (Rt) — Hq — Rf • Q . 



(3) 



Here t is time, Rt = {R\ , . . . ) denotes a set of time de- 
pendent external parameters, (the generalized "displace- 
ments") Q = (Q 1 , ■ . ■ Q ) is the set of conjugated sys- 
tem observables (the generalized "forces"), and "•" de- 
notes the scalar product. We assume that the bath is 
ideal, meaning that it is has infinite heat capacity C and, 
accordingly, cannot change its temperature T upon in- 
jection of finite amounts of energy. For convenience we 
introduce the notations 
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to denote the Gibbs equilibrium of the total system 
at fixed parameters R, and the corresponding partition 
function (Tr is the trace over the total system). We as- 
sume that at t — the total system is in the equilibrium 
£r • We next consider some system observable O, and 
ask how its expectation value at time t = r, 



(0 T ) = Tr q^Ot = Tr g r O, 
deviates from its equilibrium expectation value 
(0) e R q T =Tr^O. 



(5) 
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Here, O r denotes the Heisenberg picture O t = U^OUto, 
g t is the total system density matrix at time t: Qt = 
UtoQ^UtQ , and U T t denotes the quantum time evolution 
operator from time t to time r, generated by the total 
Hamiltonian 7f.(R t ). 

Using the cyclic property of the trace operator, and 
the property u\ t U T t = 1, one can prove the following 
nonequilibrium identity 
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d£Rt ■ Qt , (8) 



W dis = W - AF 
Eq. reads: 



dtR* ■ [Q t - (Q) c R q J , (9) 
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The operators W and PV<jis d° n °t correspond to any 
quantum- mechanical observable (22l . |30| . but approach, 
- in the classical limit -, the exclusive work, w, and dis- 
sipated work, w ~ AF, respectively [2^]. Under our as- 
sumptions that the bath has infinite heat capacity, the 
nonequilibrium expectation (0 of Wdis vanishes in the 
adiabatic limit (for a discussion of the scaling of Wdi s 
with the bath size in a classical set-up, see [32j). Since 
the expectation of Wdis is given by the Kullback-Leibler 
relative entropy between g T and g^ [33-36], this also 
means that in the adiabatic limit g T — > g^ ■ As the 
driving speed increases, the actual state g T lags more and 
more behind the "reference" equilibrium state (32l |. 

In the following we consider slow (quasi-adiabatic) 
driving, and accordingly expand Eq. (|10[) to first order 
in Wdis (here slow means that the characteristic time of 
variation of the driving is small compared to the char- 
acteristic time of relaxation to the Gibbs equilibrium). 
Following the method outlined in Ref. [27j], we use the 
operator expansion: 
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Setting A = —H T (R T ) and e = —Wdis, we arrive, to first 
order in Wdi s , at the result 

(AO T ) = (O t ) (0) C £ (12) 

= - f T dt f du(0 T e- un -^ AQ t e uU ^^ ■ R t )^ . 
Jo Jo 

Using the cyclic property of the trace operator and uni- 
tarity \j\^J T t = 1, we rewrite the integrand in Eq. (fl2|) 
as: 

TrgrOUroe-^^UrtAQuUe^^Ulo ■ R t . (13) 

Since this expression is already of first order in Wiis we 
can replace the exact density matrix g T with the approx- 
imate equilibrium density matrix g^ . The next, crucial 
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assumption is that the correlation function in Eq. f| 12[) 
decays quickly compared to the time scale of variation 
of R t , which, in fact was assumed to be very large. Un- 
der this assumption, one can approximate the exact time 
evolution operator U T t with the time evolution at frozen 
R = R T : U Tt ~ e -<«(R T )(*-T)/ft j replace R t by R T , to 
arrive at: 

(AO T ) = -E [ T ^S(t-r)K (14) 

*«5(t) = / du(AO- lhu AQl)% T . (15) 
Jo 

Note that (t) is the quantum equilibrium correlation 
function between O and AQ l (i.e., the relaxation function 
Q) calculated with respect to the equilibrium state and 
propagator at fixed R = R r . 



III. GEOMETRIC FRICTION AND 
GEOMETRIC MAGNETISM 

The theory applies regardless of the number N of driv- 
ing parameters. Geometric magnetism only appears in 
the case where there are at least N > 2 driving parame- 
ters. 

Choosing O as the i-th component of the force, Q l , Eq. 
(fT2j) becomes, using vector notation: 



(AQ T ) = -K(R T )R r 



(16) 



where K(R r ) is the N x N conductance matrix whose 
elements are the integrated force-force equilibrium corre- 
lation functions: 

K*(R T ) = [ T dt f du(AQ j _ iha AQ k t _ T r* T (17) 
Jo Jo 

= jfa MQUr m Qtr)Z ~ (Q j )Z(Q k )n T ) ■ 

(18) 

Note that the r.h.s of Eqs. (p~4| fT6|) are geometric, i.e., 
they depend on time only through the time dependent pa- 
rameters R r . Following Berry and Robbins [7| we assume 
for simplicity a parameter space of dimension N = 3 and 
rewrite Eq. (|16p in vector notation as: 
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-K S (R T )R T - B(R r ) x R T 



(19) 



The first term, stemming from the symmetric part K s 
of the conductance matrix K, is geometric friction, and 
the second term, stemming from the antisymmetric part 
K A of the conductance matrix is geometric magnetism. 
The field of geometric magnetism, B, has components 
k £ ijkK 3 ' [sijk is the Levi-Civita tensor], 
and reads, in vector notation: 



B(R T ) = -1/2 / dt f du(Q- ihu x Q t _ T )^ 
Jo Jo 



(20) 



The theory may be generalized to account for general 
nonlinear driving, i.e., to system Hamiltonians i?(R t ) 
not necessarily of the form in Eq. ([3]). In the general 
case the operator W reads [3]: W = J Q T dtR t • VH t (R t )- 
Accordingly, the general theory is obtained by replacing 
everywhere Q with V-ff . A main result, therefore, is that 
the field of geometric magnetism emerges as: 

poo t>P 

B(R) = y 2 / dt du(VH_ lhu xVH t ) c £. (21) 
Jo Jo 

where for simplicity we have changed the integration do- 
main from [0,t] to (— co, 0], we have used the fact that 
the antisymmetric part of the relaxation function is odd 
under t — > — t [|[ to express the result as an integral from 
to oo, and we have dropped the time label r = in R. 
Eq. (|2ip is the open quantum system version of Berry 
and Robbins expression (p} for the geometric magnetism 
of an isolated classical chaotic system. 

It is worthwhile to re-express Eq. (|2Tj) in terms of the 
symmetrized force autocorrelation function 



v" k {t) = y2({d k H,d j H t })%, 



(22) 



where {■, •} denotes quantum anti-commutator and dj - 
d/dRj. To this end we rewrite it in tensor notation as 



/•OO 

&(il) = (-1/2) 5> ifc / dt*£(t) ; 



(23) 



where 3>j|(£) is the relaxation function between djH 
and dkH, sec Eq. (fT5|) . According to the fluctuation- 
dissipation theorem [8j 



/oo 
dt'T(t-t^f k (t') 
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(24) 
(25) 



hence 
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dt 



jk 



or, in vector notation: 



B(R) = i/4 / dt 



dtT(t-t')*? k (t'), 

-oo 

(26) 



dt'T(t - t') 



{VH x VH t < - VH t > x VH) C £ . (27) 

Eq. (|23[) can be re-written in a remarkably simple 
form: 



Bi(R) = (-1/2) ^e yfc $f fe (0) ; 



(28) 



jk 



where the symbol $^|(s) = J °° die st( l > ||(t) denotes the 



Laplace transform of $^|(t). 
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Classical limit 



expressed in terms of the reduced density matrix 



Eq. (fT6| holds true also classically. The derivation 
can be repeated following the quantum derivation given 
above, allowing observablcs to commute. As a result, 
the quantum thermal correlation functions have to be 
replaced by the classical expressions Q , so that the clas- 
sical geometric magnetism reads: 



S cl (R) = P/2 / dt(VH x VH t )% 



(29) 



This result may also be obtained by taking the limit H — > 
of Eq. (f2Tj) . Alternatively one can take the limit H — > 
of Eq. 423). In this limit T(t) -> (36 {t) @, where S 
denotes Dirac's delta function, and observables commute 
V.HV x VH -> -VH x VH V . 



IV. OPEN SYSTEM DYNAMICS 

The geometric magnetism (and the geometric friction 
too) may be recast in the more familiar language of dis- 
sipative open system dynamics [37l |38| , in terms of the 
system reduced density matrix 



(30) 



where Tr# denotes the trace over the bath Hilbert space. 
The linear response of the force, which defines the con- 
ductance matrix K (hence the geometric friction and the 
geometric magnetism), may be written as: 



(AQ T ) = Trp T Q - Trp^Q 
= Tr sP s T Q-(Q)Z S : 



(31) 



where Trg denotes the trace over the system-5 Hilbert 
space and denotes ex pec tation over the equilib- 

rium reduced density matrix |39j 



cq,S 
PR 



Tr B e 
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,-0H*(R] 



/Zr , (32) 



where H*(R) is the Hamiltonian of mean force and 
Zr = Ttse ^ is the partition function of an open 
quantum system [2^, I37l - l39j . In the case of weak cou- 
pling the Hamiltonian of mean force reduces to the sys- 
tem Hamiltonian Hg. 

The element K^ k of the conductance matrix may be 
experimentally/numerically obtained by driving the sys- 
tem with a small constant velocity in the j direction, 
Vjj and measuring/computing the reaction force in the 
k direction: 



K> k (R T ) = - Tr s tf ,j Qk - (Q, 



\eq,S 
>R T 



/Vj 



(33) 



where we have introduced the notation J to denote 
the reduced density matrix resulting from the perturba- 
tion Vjj . Accordingly, the geometric magnetism may be 



TrsP^Qk-iQkTjT IV, ■ (34) 



eq,£> 
R 



As illustrated above the geometric magnetism may be 
accessed also by calculating the equilibrium force au- 
tocorrelation function. Although quantum correlation 
functions can not in general be expressed as expectations 
over the reduced density matrix, they are "open quantum 
systems" objects that depend explicitly on bath proper- 
ties, notably the bath spectral density, see Eq. (|39[) be- 
low. In particular, it should be noted that exact open 
quantum system dynamics generally is (i) neither lin- 
ear [4(| (ii) nor can it be described by trace preserving 
completely positive maps [U |42[. Attempts to resort 
to approximations (e.g. Markov, and rotating wave) to 
express the correlation functions in terms of Markovian 
dynamics for the system observables, may lead to results 
which contradict basic principles such as the fluctuation 
dissipation theorem [43j , and therefore to non- negligible 
errors in the evaluation of geometric friction and mag- 
netism. See Ref. [44[ for a recent example of the drastic 
effects that even good approximations may have on the 
calculation of geometric phases. Therefore a very special 
care must be paid when employing such approximations 
in this context. 



V. ILLUSTRATION: THE DAMPED CHARGED 
HARMONIC OSCILLATOR IN A MAGNETIC 
FIELD 

As an illustration of the theory we consider a quan- 
tum damped charged harmonic oscillator of mass m and 
charge q transported along a path R t in presence of a con- 
stant magnetic field B. Adopting the Caldeira-Leggett 
model of quantum Brownian motion (45j . the system, 
bath and coupling Hamiltonian read: 

H(R t ) =(p- qA) 2 /(2m) + muj 2 x 2 /2 - muj 2 x ■ R t 

N 

H B = Y.^ p n/ m n + m n u>l$l]/2 



N 



N 



H S b = -x ■ Cn£ n + x 2 ^2 c„/{2m n ujl) . (35) 



Here x,p, to denote the harmonic oscillator position, mo- 
mentum and frequency, respectively. £ n , P n , m n , co n de- 
note the nth bath's oscillator position, momentum, mass 
and frequency, respectively. The symbol c„ denotes the 
linear coupling constant between the harmonic oscillator 
and the nth bath's oscillator. The symbol A denotes 
the vector potential. Note that according to Eq. © 
Q = muj 2 x. Assuming an initial global Gibbs distribu- 
tion and adopting the Fcynmann- Vernon path integral 
approach (37| one arrives, after integrating out the bath's 
degrees of freedom, at the following generalized quantum 
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Langevin equation for the charged oscillator's position 


mx t + / dt'rj(t — t f )xf — qx t x B + mcu 2 x t = F t + ft-> 

J — oo 

(36) 

where 77(f) is the friction kernel, f t — muj 2 R t is the 
externally applied force, and F t is the stochastic force. 
Without loss of generality we shall assume that B points 
in the z direction, B = Bz. Since the motion in the 
z direction is decoupled from the motion in the x and 
y directions, the xz and yz relaxation functions vanish, 
implying that the geometric magnetism is also directed 
in the z direction. Further, due to spatial homogeneity, 
geometric magnetism does not depend explicitly on the 
position R t . That is, B = Bz. From the compact ex- 
pression the strength of the geometric magnetism 
reads 



B 



1/2^(0) -**„«))] = -^(0) (37) 



where we have introduced the notation , for the anti- 

symmetric component of <fr xy and $> xy (s) f° r ^s Laplace 
transform. Following [46| the Laplace transform of the 
antisymmetric part of the response function reads 



&Sv(«) = 



(mco 2 ) 2 qBs 



sfj(s)] 2 + q 2 B 2 s 2 



(38) 



where [■, ■] denotes quantum commutator, and fj(s) is the 
Laplace transform of the bath friction kernel. Its form 
depends on the bath spectral density. For an ohmic bath 
fj(s) is constant. As compared to Eq. (2.15) of Ref. [46| . 
we have in our Eq. (|38[) an extra factor (mto 2 ) 2 stemming 
from our definition of a xy in terms of Q x = mu) 2 x, Q y = 
mu) 2 y, rather than x,y. 

Since $> xy = f°° dt'a xy (t') [1] we have 



5" (0) 



(mu 2 ) 2 qB 



sf)(s)] 2 + q 2 B 2 s 2 



{mu 2 ) 2 qB 



\mu> 



s 2 + st)(s)] 2 + q 2 B 2 s 2 ' 



(39) 



where a% y (Q) = due to the fact that at equal times Q x 
and Q y commute. Regardless of the bath spectral den- 
sity, the friction kernel 77(f) vanishes at large times, hence, 
according to the final value theorem lim s _>o s v( s ) = 0- 
Using Eq. ([37)1 one finally obtains the result 



B = qB, 



(40) 



which evidently holds both classically and quantum- 
mechanically. Apart from the charge q, geometric mag- 
netism is nothing but the physical magnetic field in this 
case. The factor q stems from the fact that the geomet- 
ric Lorentz force in Eq. (|T9"|) reads —BxR, whereas the 
Lorentz force reads — qB x v [where v = (x, y, z)]. This 



very same result was found also in Refs. [43, |48[ for the 
case of an isolated classical or quantum harmonic oscil- 
lator. Our result (|40|) conveys the non-trivial knowledge 
that this continues to hold also for an open classical or 
quantum harmonic oscillator. That is, the presence of a 
bath does not destroy the geometric magnetism, in fact 
it does not minimally alter it in this case. Analogous cal- 
culations involving the symmetric part of the relaxation 
function lead to the result that the geometric friction is 
given by the time integral of the friction kernel J Q dtrj(t). 

It is noteworthy that the case of geometric magnetism 
is distinct from the case of standard equilibrium diamag- 
netism, which is absent in the classical limit of open sys- 
tems and reveals itself at the quantum level only, see 
Bohr- van Leeuwen theorem W3, Mm . 



VI. CONCLUDING REMARKS 

We have derived a general expression for the field of ge- 
ometric magnetism in open quantum systems, Eq. (|2ip . 
possibly coupled strongly to the environment. This gen- 
eralizes the expression (J]) of Berry and Robbins [D, 0] 
which refers to closed chaotic classical systems. It is 
worth noticing that, contrary to the case studied by Berry 
and Robbins, here no assumption of chaotic dynamics of 
the driven system H(R t ), which may well be integrable, 
is made. It is the presence of the thermal bath Hb and 
the coupling to it, Hbsi that provide the necessary de- 
gree of chaos for the development of a response theory a 
la Kubo. It is however important to remark the differ- 
ences between the presently developed theory and that 
of Kubo. This is best seen by confronting Eq. (fl"4"|) with 
Kubo's formula 



(Or) - (0)2, = / dt / dw(AO_ ift „AQ t _ T )^ ■ R t . 
Jo Jo 

(41) 

Note the prominent difference that Kubo's formula (|4"Tj) 
gives an expression (linear in R) for the difference be- 
tween the nonequilibrium expectation of O at time r, 
and its equilibrium expectation at time 0, whereas the 
present formula (|14[) gives an expression (linear in R) for 
the difference between the nonequilibrium expectation of 
O at time r, and its equilibrium expectation at the same 
time t. Thus in Kubo's theory the small parameter is the 
strength of the driving whereas in our theory the small 
parameter is the speed. Both formulae (fT4l l4Tj) yield the 
response in terms of equilibrium correlation functions. 
While Kubo's formula involves the correlation between 
O and Q (the response function), our formula involves 
the correlation between O and Q (the relaxation func- 
tion). Note that Kubo's formula (|4"Tj) follows from an 
exact fluctuation relation 



(O r e 



TO , 



(42) 



that looks very similar to our starting Eq. ((?]) (2(| [27| ■ 
The differences are that (i) the r.h.s. is evaluated at Rq 



6 



in Eq. (|42|) . while it is calculated at R T in Eq. (J7J, 
(ii) Eq. does not involve the free energy difference 
AF, which instead appears in Eq. (0, (hi) in Eq. (j42|) 
the unperturbed system Hamiltonian Hq appears instead 
of the total Hamiltonian TL(Rt) appearing in Eq. (J7J. 
These complementary expressions (0 and (|42|) are cus- 
tomarily referred to as "inclusive viewpoint" and "exclu- 
sive viewpoint" fluctuation relations, respectively. Inter- 
ested readers can find accounts of the importance of these 
viewpoints in the theory of nonequilibrium fluctuations 
in Refs. [H |2j|. Just like Eq. (JUJ) allows one to 
obtain Kubo's formula (|4"Tj) and the whole hierarchy of 
higher order nonlinear responses, so does Eq. ([TJ allow 
to obtain the adiabatic linear response relation (fRij) . as 
well as the higher order terms in the adiabatic expansion. 
An interesting open question is whether and under which 
conditions geometric forces appear in those higher order 
terms. 

Our main result, Eq. (|21j) . provides a straightforward 
way to define the Berry phase of an open quantum sys- 
tem. Just like the surface integral of the classical two- 
form ([T]) provides a generalization of Berry phase for 
chaotic classical systems [1, Q , so does the surface inte- 
gral of the geometric magnetism (|2T|) provide an analogue 
of the Berry phase of open quantum systems, reading 

7 = J B • dS . (43) 

This so defined phase 7 would in general differ from those, 
equally sound and useful, expressions of a Berry phase in- 
troduced for open systems in the prior literature [TlTl2ll ] . 
In full analogy with the original Berry phase, 7 is geo- 
metric, that is it depends only on the path described by 
the driving parameters. It vanishes for a path enclosing 
no area, and it vanishes in the case when the system dy- 
namics are time-reversal invariant, i.e., when for any t, 
QH(Rt) = "H(Rt)<d. Here O is the anti-unitary time- 
reversal operator which reverses momenta and keeps the 
spatial coordinates and all external parameters (possibly 



including physical magnetic fields) unaltered 111. 12211 This 
is so because, due to Onsager-Casimir relations [51j, the 
conductance matrix K would be symmetric in this case, 
hence, the geometric magnetism 3 would vanish. 

Our simple example of a quantum harmonic oscillator 
transported along a path already shows that the pres- 
ence of an environment does not destroy geometric mag- 
netism. In fact, in this specific (linear) case the geo- 
metric magnetism is given by the actual physical mag- 
netic field, exactly like in the isolated case [13, ED]. For 
nonlinear systems the difference between the real and 
geometric magnetic fields could be detected, as well as 
the difference between quantum and classical regimes. 
However, the quantum-mechanical treatment of nonlin- 
ear open systems constitutes an ambitious challenge be- 
cause in this case the system evolution cannot be handled 
analytically in an exact manner. This challenge, in prin- 
ciple, could be approached numerically, for example, (i) 
by resorting to the Floquet-Markov formalism [52j , under 
the assumption of weak system-bath coupling; or (ii) by 
following the Feynmann- Vernon path integral formalism 
[29l Wf\ , to calculate the reduced density matrix numeri- 
cally, through stochastic unraveling of the corresponding 
influence functional [53j . 

Geometric magnetism is at the basis of a currently 
growing experimental activity aimed at producing artifi- 
cial gauge fields in thermally isolated cold atomic gases 
|54T - [58| . The present theory opens the possibility of en- 
gineering synthetic gauge fields also in the presence of a 
thermal environment, via our general expression (1211) . 
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